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1 . Introduc t ion 

Carrier Sense Multiple Access with Collision Detection (CSMA/CD) is a 
channel access protocol for packet broadcasting local area networks. This 
protocol allows user stations to contend for time on a shared channel by 
multiaccessing it in a random fashion. In an attempt to reduce the frequency 
with which broadcast message packets collide on the shared channel, Ethernet 
(Metcalfe and Boggs 1976) and IEEE 802 (IEEE 1985) use a technique known as 1- 
persistent CSMA/CD with binary exponential backoff (Goodman et al . 1988). An 
alternative technique proposed for reducing the frequency of collisions is p- 
persistent CSMA/CD (see, for example, Stallings 1987 and Takagi and Kleinrock 
1985a); this protocol is considered in this paper. 

The dynamic behavior and performance of unbuffered CSMA/CD systems has 
been studied extensively (e.g., Lam 1980, Takagi and Kleinrock 1985a, and 
Tobagi 1980) . These studies assumed that there are an infinite number of 
network users and each user can have at most one packet in the message queue 
at any time (called an infinite source model). This simplifying assumption has 
been relaxed in recent studies on CSMA/CD systems: Apos tolopoulos and 
Protonotarios (1986), Goodman et al. (1988), Hastad, Leighton and Rogoff 
(1987), Park and Bartoszynski (1990a), Park and Bartoszynski (1990b), Takagi 
and Kleinrock (1985b), and Tasaka (1986). In these studies, a finite number of 
users are considered to have message queues with finite or infinite capacity 
(called a finite source model or a buffered CSMA/CD system) . 

Among these studies on buffered CSMA/CD, the following authors considered 
p-persistent CSMA/CD. Under the assumption of finite capacity queues, 

Apos tolopoulos and Protonotarios (1986) modeled the message queueing process 
as a two-dimensional semi-Markov chain. With the state space reduced to a two- 
dimensional space, it was necessary to employ an iterative approximation 
procedure to obtain the mean packet delay in steady state. The transmission 
probability p was arbitrarily set equal to the inverse of the number of busy 
users in the beginning of each contention period. Thus the dependency of the 
mean packet delay on the control parameter p was not investigated. Takagi and 
Kleinrock (1985b) developed a stationary joint probability generating function 
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(PGF hereinafter) of queue lengths for a system with two users, and obtained 
an explicit formula for the mean packet delay based on that PGF. An exact 
joint PGF of the queue length vector for a system with more than two users was 

t 

first developed by Park and Bartoszynski (1990b). Further, they obtained 
stability conditions analytically based on the PGF. 

In implementing p- persistent CSMA/CD the major problem is to determine 
the value of p so as to minimize the response time for users while maintaining 
stability of the system. In this paper an approximation procedure is developed 
to compute the optimal p which leads to a stationary operating mode and at the 
same time yields the minimum mean packet delay in equilibrium. An optimization 
of the p-value with respect to the response time has not been studied for a p- 
persistent CSMA/CD system with more than two users , though it has been carried 
out for slotted ALOHA by Saadawi and Ephremides (1981) and Sidi and Segall 
(1983) . 

In the next section, the p-persistent CSMA/CD channel access protocol is 
described and the queueing dynamics under this protocol is modeled as a multi- 
dimensional semi-Markov chain. A stationary PGF of packet backlog, i.e., the 
sum of all queue lengths, is obtained in Section 3 based on the joint PGF of 
the queue length vector. In Section 4 the Kolmogrov forward equation is 
derived applying an M/G/l approximation to the stationary PGF of packet 
backlog, and an iterative procedure is developed to evaluate the mean packet 
delay in equilibrium and optimize p with respect to the mean delay. A numeric 
analysis is conducted in Section 5 and the analytic results are compared with 
simulation results demonstrating the accuracy of the proposed approximation 
scheme. Section 6 contains concluding remarks. 

2. A Semi -Markovian Model of the Message Queueing Process 

A brief description of the p-persistent CSMA/CD channel access protocol, 
which is sufficient for the reader to follow the subsequent analysis, follows. 

The system consists of a single server (the channel) and multiple, 
inf ini te - capacity queues of customers (message packets). The number of queues, 
denoted by m, can be as large as 200 in local area networks. The channel time 
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is slotted with the slot size being the maximum propagation delay, and a slot 
size is chosen as the unit of time. Discrete time index t-0 , 1 , . . . is used to 
denote slot boundaries. 

If the channel is busy no users attempt transmission (due to carrier 
sensing ) until ongoing transmissions are completed and the channel becomes 
empty again. Once the channel becomes empty, every busy user (i.e., a user 
with a nonempty queue) persistently attempts at every t to transmit a message 
packet from the top of its queue. In a transmission attempt a user samples a 
random number from the uniform distribution over [0,1]. If the number is 
smaller than p, it starts transmission of a packet; otherwise, it suppresses 
the start of transmiss ion . Thus whenever the channel is empty all busy users 
contend to seize the channel with an equal chance. If only one among all busy 
users samples a random number smaller than p at time t, the user starts 
transmitting a packet at t and the packet is successfully transmitted. If two 
or more users start transmitting packets at the same t, packets collide, the 
collision is detected, and all collided transmissions are aborted before 
transmissions are completed. 

It is assumed that in each slot at most one packet arrives at a queue 
(which is reminiscent of Poisson process assumption) . The probability that a 
packet arrives at a queue during a slot is denoted by A. The numbers of 
packets arriving in different slots at the same user are independent, and the 
numbers of packet arrivals at different users in the same or different slots 
are independent. This message input process (called a geometric arrival 
process ) implies that queueing processes at different users are symmetric; 
viz., queue lengths at individual users are exchangeable random variables 
(Billingsley 1979) . 

Following Takagi and Kleinrock (1985a), the channel state dynamics are 
modeled as depicted in Figure 1 (note that users 3 to m - 1 in the figure are 
assumed to be idle throughout the portion of time shown in the figure) . The 
channel state alternates between idle periods and busy periods . During an idle 
period all users have empty queues, and during a busy period there is at least 
one busy user. With the assumed message input process, durations of idle 
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periods are independently and geometrically distributed; hence the beginning 
of an idle period is a system regenerative point. A busy period is divided 
into a number of sub (busy ) periods , each in turn consisting of a transmission 
delay followed by a packet transmission period. 

Figure 1 about here 



During a transmission delay no packets are transmitted, even though the 
channel is empty, since all busy users continually sample random numbers 
exceeding p. A transmission delay terminates at time t when one or more busy 
users sample random numbers smaller than p and start packet transmiss ions . The 
duration of a successful packet transmission is proportional to the length of 
the packet transmitted. Packet lengths may vary. However, for analytic 
convenience, we use the expected value of packet lengths and denote it by Z. 
The duration of an unsuccessful transmission, which is the time to detect and 
resolve a collision, is represented by Z' . In a baseband system like Ethernet, 
it takes at most two slots to detect a collision. Thus Z' is usually much 
shorter than Z. A p- persistent CSMA/CD system is completely described by the 
set of parameters, { m , p , A , Z , Z ' ) , at least for the purpose of our analysis. 

Throughout the paper we use the following convention in defining 
notations: A boldface capital letter indicates an m-vector or a function with 
vector arguments, except that E and P stand for expectation and probability, 
respectively. An italic capital letter with an underlined superscript, i, 
indicates the ith element of an m-vector, and an underlined italic capital 
letter signifies the sum of all elements in an m-vector. An italic capital 
letter without a superscript or an underscore signifies a set or a function 
with scalar arguments . 

Denoting the set of network users by 1= ( 1 , 2 , . . , m) , we define Q^=* [Q“, 
i£l] to be the vector of queue lengths at c, and ^=1^ to P ac ^ eC: 

backlog at t. The process ( Q ^ ) is characterized as an irreducible, aperiodic 
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discrete- time Markov chain whose state space is the m- dimensional vector space 
O-Uwj, iel]: uk> 0 , integer}. 

A semi-Markov chain model of the queueing process is constructed as 
follows (which is more accurate than the semi-Markov chain model presented in 
Takagi and Kleinrock 1985b): The beginning and the end of each idle period are 
regarded as embedded Markov epochs. Within a busy period, the boundaries of 
subperiods are chosen as embedded Markov epochs. Denoting these epochs by 
k-0,1,.., with k-0 at t-0 , we have an embedded Markov chain {Q^} which is a 
multidimensional random walk. In Figure 1, the embedded epochs are shown by •, 
and the packet backlog at those epochs, are illustrated. 

3. Stationary Probability Generating Function of Packet Backlog 

The major difficulty in analyzing the queueing process (Q^) lies in that 
queues at different users interfere with each other through the shared 
channel. This problem of interacting multiple queues cannot be solved using 
classical queueing theory. There are no analytic results in the literature for 
evaluating the mean packet delay in such a multidimensional queueing system, 
and even approximations are difficult to obtain when m> 2 (see Kleinrock and 
Yemini 1980) . In this paper we base the derivation of the mean packet delay in 
equilibrium on the stationary PGF of IQ^) • By the aggregation of (Q^J into a 
one-dimensional process (2^) we . of course, lose information about behaviors 
of individual queues. The approximation needed to overcome this loss of 
information can, however, be localized in the delay analysis (as will be shown 
in Section 4) so as to generate accurate results (verified by simulation in 
Section 5) . 

To describe the process (2^1 we first define the following random 
variables. Notations for each element in a random vector and the sum of all 
elements in a random vector are defined according to the convention mentioned 
earlier . 

X(n) : numbers of packet arrivals at respective queues over an n-slot period. 
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Z k : numbers of packets that are successfully transmitted from respective 

queues at the end of the subperiod which begins at the kth epoch; Note 
that either Z^=0 (O' vector ) or *1 for some i and Zr *0 for all j>i. 
R k : duration (in number of slots) of the transmission delay in the sub- 
period beginning at the kth epoch. 

L k : duration (in number of slots) of the packet transmission period in the 

subperiod beginning at the kth epoch. 

U k Q k+ 1" Q /c' 

B k { iel : 2jr >1), i.e., the set of busy users at the kth epoch. 

S^(i?k) {iel: 2p=0 , X-(i?^)>l } , i.e., the set of users who are idle at the kth 

epoch, but become busy by the end of the transmission delay in the 
subperiod beginning at the kth epoch. 

The process (2^) can now be recursively defined as 

2k+ 1“ where 0; X(R k +L k ) -Z k otherwise. (1) 

In this equation X(l) represents random variable Z(l) with the constraint 
X(l)> 1. Under the geometric input process, X(l) has a binomial distribution 
with parameters m and A. Therefore, if 2^=0 (i.e., if the kth epoch is the 
beginning of an idle period) , 2^ + ]_ follows a truncated binomial distribution 
given by 

?(Q k+l =q\Q k =0) - 1< q< m. (2) 

Note that Q k +^ gi ven 2^=0 is the packet backlog at the start of a busy period. 

Next, if Q k > 1, the random drift U_ k is the sum of random variables X(R k ) , 
X(L k ) and -Z^. Unfortunately, the latter variables are stochastically 
dependent on each other in a complex way. First, L k (and hence Z(L^)) depends 
on Z^. Namely, L k =2 9 if Z^=0 and otherwise (i.e., if Z^-l) . Random 

variable Z^ in turn depends on \& k \ and since the success of a packet 

transmission depends on the number of busy users at the end of the transmis- 
sion delay, which is the sum of \B,\ and \B',(R, )|. The latter, \Bl(R, )\, 
depends on X(R^) ■ Further, R^ and are jointly dependent on |3^|, and 
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finally |B^| is uniquely determined given Q^, but stochastically dependent on 
2^. These stochastic dependencies are summarized in Lemma 1 below. A detailed 
analysis of these stochastic dependencies and the resultant conditional PGF of 

i 

^k+1 &i ven Qfc-1 can be f° und in Park and Bartoszynski (1990b). A derivation of 
the latter PGF and the proof of Lemma 1 are summarized in Appendix I at the 
end of this paper. 



Lemma 1 : For any k with 2^0. 

/ 

E(t/ k |2 k )- EE|x(i? k )+EE 



E[i(l k )|z fc ]-Z, 



IWI 



R k .l(R k ) 



l«. 



An optimization of the p-value with respect to the mean packet delay has 
been studied for slotted ALOHA, as mentioned earlier in Section 1. Compared 
with slotted ALOHA, the delay analysis and optimization of the p-value is 
significantly more complex for p- persistent CSMA/CD. The source of additional 
complexity is the stochastic dependencies among the random variables as 
described in Lemma 1. These stochastic dependencies stem from the carrier 
sensing and collision detection mechanisms which have been added to the 
slotted ALOHA protocol in devising p- persistent CSMA/CD (see Szpankowski 1988 
and Tsybakov and Mikhailov 1980 for a Markov chain model of slotted ALOHA with 
buffered users). In this paper the stochastic dependencies among L^, 

X(k^) and | £^| are exactly incorporated in the delay analysis. An 
approximation is introduced to simplify the dependency of |£^| on 2^ * 

The process (2^) is a random walk with random drifts (2^)* The process is 
uniformly downward bounded meaning that P (U^<-1)^0 for all k. It also 
possesses an important property called the bounded homogeneity , which is 
described in the following Lemma (see Appendix I for the proof): 



Lemma 2 : For any k, k' , 2^ and 2^/ such that | BjJ ~ | | , 12^)“ 

P(^ k ,=u |2 k ,). 

Define F k (s) to be the PGF of Q k - In view of (1) we write 
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F k+l (s >- F k+ i(s|2 t - 0 > p <2k- 0 ) * F i+ i <s l2 k >0)P(2 t >0). 



(3) 



LetCing F k+ i ( s | • we obtain from Equation (2) 

H(s)- { (As+1-A) n, -(1-A) m )/(1-(1-A) ra ) . 
By Lemmata 1 and 2 we may write 



(4) 



f 'k + l (s l2 k >0)- E 



2|t„J 

s EEs s 



IS. I 



2 t >0 



^ «o min (q,m) ® 



^ s U P(£/ k =u | |B k |=b; P(|B k |=b|2 k =q) PCS*-?). (5) 



P(2 k >0) q-1 5=1 u=-l 

Define C^(s \ \ B ^ | =5) := _-^s u P (I/^=u | | | =b) , i.e., the conditional PGF of 

given By Lemma 2, for any b* 0, G(s|b):= G^(s | | | =*fc>) is independent 

of k and . By setting s— =s for all iGl in Equation (A4) in Appendix I, we 
get 

m - hi 

\s 



m- b 

G(s|b)= ^ 

J- o 



r 

.-a 

i 

e 




1 j , 





( 1 ■ P ) (As+l-A) - (l-A)J 



J 



X 



5* 1 (As+l-A) ra ' e (b+j)p(l-p) b+J '- 1 +( As+l-A)^' {l-(l-p) b V.(b+j ) p( 1 -p) b+ J-l ) 



X 



n-0 






1- (l-p) b+n (l-A) m * b ' n (As+l-A) b+n j 



-11 



( 6 ) 



Note that index j in Equation (6) represents the number |B^(£^)|. 

By substituting (6) into (5), (4) -(5) into (3), and then by passing k to 
oo in (3), the stationary PGF of (2^1, denoted by F(s), is written as 
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min(g , m) 



(7) 



F(s)~ tf(s)P(fi-0) + ^ s q 

<7-1 



b-1 



C(s|b)P(|B|-b|2-q)P(2-q). 



As in (7) we simply drop the time index k from an argument to indicate its 
quantity in steady state. 



4. Delay Analysis 

In the sequel, we write F^(c) to denote the jth partial derivative of 
F(s) with respect to s, evaluated at s=c . We need to evaluate F^^(l) to 
obtain the mean packet backlog at an embedded epoch in steady state, which in 
turn is used to evaluate the mean packet delay. Analytic evaluations of F(s) 
and its derivatives are difficult to obtain mainly because the conditional 
probability P ( | B | =fc> 1 0=q) cannot be evaluated explicitly in the presence of 
interactions among queues. Therefore, we introduce an M/G/l approximation in 
formulating this probability. 

In this approximation, we assume that in steady state the message queue 
at a user behaves like a discrete M/G/l queue independently of other queues. 
The queue at a user in the CSMA/CD system indeed has a single server (the 
channel) and geometric arrivals (a discrete version of Poisson arrivals) of 
message packets. But how about the service time? The service time for a packet 
begins at the point when the packet is promoted to the top of its queue. The 
service time is defined, in the context of an M/G/l system, to consist of: (i) 
a number of sub (busy )periods , since that beginning point, during each of which 
the packet is not transmitted (due to sampling a random number greater than p) 
or is involved in a collision; plus (ii) the subperiod in which the packet is 
finally transmitted with success. This service time is actually affected by 
other busy users' contention for the server. Furthermore, the set of busy 
users may change during the service time of a particular packet. The 
assumption of independence (or "no interference") among the queues in steady 
state is thus introduced so that service times are i.i.d. in steady state. 

With this assumption, we now may write 
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P< IB I -fa)- ^8 b a-8) m ' b 



( 8 ) 



where 9 is defined to be the steady-state probability of a user having a non- 
empty queue at an embedded epoch. 

The no interference assumption for the steady-state behavior of multiple 
queues should lead to a closer approximation when it is applied to p- 
persistent CSMA/CD than when applied to slotted ALOHA. This is because the 
former suffers less from interference (viz., collisions) than the latter due 
to additional controls, carrier sensing and collision detection . The 
approximation in (8) is motivated by the results in Sidi and Segall (1983) 
where an M/M/1 approximation is applied to the slotted ALOHA demonstrating 
good accuracy over a wide range of system parameters. It should be noted that 
the approximation is limited only to computing 9 , and the rest of delay 
analysis exactly reflects the dynamic behavior of interacting multiple queues. 

Incorporating (8) into (7) and using the equation F( 0)- P(2~0), we obtain 

F(s|0)- H(s)F(0\ 9) + V(s|0){F(s|*)-F(O|0)}, (9) 



where 



V(s 0) 



C(s 



b=l 






9 u (l-9 ) 



m- b 



( 10 ) 



Taking the first derivatives of both sides of (9) yields 



•a;, 



(s\9)- H (1) (s)F(0\9) + V (1) (s\9){F(s\9)-F(0\d) )] /{ 1-V(s | 0) ) . (11) 



By evaluating Equation (11) at 1 and using the condition F(l|0)~ l/(l|#)- 1, 
the equilibrium differential equation (i.e., Kolmogrov forward equation in 
equilibrium: see Kleinrock 1975) can be derived as 



0 = ( 1)B(0 | <9 ) + ( 1 1 0 ) (1-F(0| $) ) , 



( 12 ) 
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which can be rewritten as 



F(0\6)- V (1) (l\9)/(V (1) (l\6)-H (1) a)) . (13) 

Here and V^\l\d) are the expected drift of packet backlog during an 

idle period and that during a sub (busy)period , respectively. From (4) we get 
( 1) - mX/( 1— (l\d) is obtained simply by replacing G(s|b) by 
(l\b) in the RHS of Equation (10). (l\b) is the expected drift of 

packet backlog over a subperiod beginning with b busy users, and is derived in 
Appendix II using Equation (6) . 

In light of (12), the condition (l\d)< 0 must be satisfied for the 

system to achieve equilibrium. In Equation (10) we see that V^\l\ B) is the 
marginal expectation of G^^(l|b) over 5*1,.., m f which in turn is equivalent 
to E(£/^| £^>0) . Thus the equilibrium equation (12) implies that a necessary 
condition for stability of p- persistent CSMA/CD (i.e., ergodicity of { Q ^ ) ) is 
approximately given by E (£/^ 1 2^>0 )< 0. Furthermore, by Pakes' (1969) lemma, the 
latter condition is also sufficient for ergodicity of the one - dimens ional 
random walk (2^). it can be easily proved that { Q _ ) is ergodic if and only if 
(2^1 is ergodic. Therefore, the condition E(2^ 1 2; c > 0) < 0 can be regarded as a 
sufficient and necessary condition for stability. However, a theoretical proof 
of this assertion is difficult and has not been established in the literature. 

t 

The reader is referred to Park and Bartoszynski (1990b) for a theoretical 
analysis of stability conditions for p- persistent CSMA/CD. Some of the main 
results in that paper are summarized in the following for later reference in 
Section 5. The conjecture below has been numerically established. 

Lemma 3 : { Q c ) is recurrent if G^^ ( 1 1 b)=E(U^ | | | *5) < 0 for all 5-1,.., m, and 
only if mA[^+(l-p) ra /{l-(l-p) ra ) ] < 1. 

Conjecture : Given system parameters im,\,2,£'}, for any b there exist p, (b) 
and Py(b) such that (?^(l|b)< 0 if and only if p^(b)< p< p^,(b). 
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Further , t Q c ) is recurrent if p L (l)< p< (m) since p L (b)< p (fc') and 
p u (b)< py(b') for any b> b'. 



Recall that 6 is the probability of a user being busy in steady state, 
and that F(O|0) is the conditional probability of all users being idle in 
steady state given the value of 9. It follows that, under the assumption of 
independent queues, we must have 1-0- F(0\ 9)^ // ' m . We thus compute the value of 
9 iteratively using Equation (13) until the following convergence criterion is 
me t : 

\6 - {l-F(O|0) 1/m ) | < £ . (14) 



The function f(8):- 6+F(0 \ 6 ) l i s monotone increasing in 8 and changes sign 
once in the range (0,1) of 9 . Therefore, a bisecting search method can quickly 
find the value of 9 satisfying (14) even for a extremely small value of e. 

Let d be the value of 9 satisfying (14). We can now evaluate F^^(l|5) by 
applying L'Hospital's rule to (11) and then replacing the RHS of Equation (13) 
for F(0|$) : 



F fl ^(l|5)- \h ( 1) (1) V (2) (l\d) - H (2) ( 1) V (1) (l\S) 



2 V (1) (l\S) V (1) (l\$) - H (1) ( l)j 



(15) 



In this equation, H^ 2 ^ { 1)= m(n?-l) A^/{ 1-(1-A) ra ) from (4), and V^^(l|0) is 
obtained simply by replacing (7(s|b) by G v ^(l|b) in the RHS of Equation (10). 
A derivation of G^ 2 ^(l\b) is provided in Appendix II. 

Recall that F^ x ^( l|£) is the mean packet backlog observed only at the 
embedded Markov epoch. In order to obtain the steady- state mean packet delay 
observed at an arbitrary slot boundary t, we need to evaluate the mean packet 
backlog at an arbitrary t in steady state (see Takagi and Kleinrock 1985b for 
a similar treatment). If c is in an idle period, the backlog is zero. The 
probability that an arbitrary t is found in a busy period when the system is 
in steady state is given by (1-F(0|$)). Letting Q denote the mean packet 
backlog at an arbitrary c in steady state, we have 
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Qr ( 1 -F(0 | $ ) )S/T 



(16) 



where 7 is the expected duration of a sub (busy )period in steady state, and 5 
is the expected backlog accumulation over a subperiod in steady state. 

First, to evaluate 7, we note that the duration of a subperiod is i.i.d. 
for all subperiods starting with the same number of busy users, due to the 
bounded homogeneity (see Lemma 2). Letting 5(b) denote the expected duration 
of a subperiod starting with b busy users, we may write 7- Z^_-^5 (b)P( |B | -b) . 
Furthermore, due to regeneracy of the queueing process the duration of a 
subperiod starting with b (>0) busy users converges a.s. to E(i^+L^ | | B-^ | -b) 
provided that the process is ergodic. Consequently, for any q> 0 we have 



r 


‘ r | -i 


* 


5(b)- E(J? 1 +L 1 | |S 1 |-b)- EjF 1 +E 


E l L ik 


R 1 ,\B[(R l )\ 



l^l-b 



(17) 



reflecting the stochastic dependencies among R^ t and Z ^ (see Lemma 1). 

Define T(s|b) to be the conditional PGF of (R^+L^) given |B^|-b>0. This 
PGF and 5(b)= T^^(l|b) are derived based on (17) in Appendix III. Using (A9) 
in the Appendix and Equation (8) we get 



7= 



m 

^ r^dlb) 

6-1 



(18) 



Turning to S in (16) , let us first define 0(q) as the expected backlog 
accumulation over a subperiod beginning with q outstanding packets so that 5“ 
2^_i^(q) P(2“^7) ♦ To derive /3(q) for q>0 , we define £(y|q f b) to be the total 
time (in number of slots) that (2 t ) spends in state y during the subperiod 
which starts with q packets waiting at b (>1) busy users. Due to bounded 
homogeneity , f(y|q,b) is i.i.d. for all subperiods that start with q waiting 
packets and b busy users. Further, it converges a.s. to its expected value if 
the queueing process is ergodic. Therefore, for any q> 0 we have 
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® rain (q,m) 

?(<})- ) l Y E?(y|q.b) P(|S|-b|2-q). (19) 

y-q 5*1 



The process (2 C ) is non- decreas ing during a subperiod and the packet arrival 
process is time homogeneous. Consequently, the average backlog at an arbitrary 
slot boundary in the subperiod which starts with q waiting packets and 5 busy 
users, or 2^ -c? yE^ (y | q , 5)/<5 (5) , is estimated by q+( . 5 )m\6 ( 5) . We thus can write 

<x> oo min(q,m) 

S- \ 0(.q)?(Qrq)- ) ) 6(5) (<?+(• 5)mA5(5)) P( | 5 | -5 | 2=q) P(£-q) . (20) 

q»l 5=1 



From (9) we see that S®, qP ("2=q)= d[ V(s | 0 ) {F(s I 0 ) -F(0 I 0) ) ]/ds j ,= F^' ) (lU) + 

q-1 | s*l 

V ^ ^ ( 1 | 0 ) { 1 -F(0 | 9 ) } . Incorporating the latter and (8) into (20) we get 



m 

S= T F^^(l|$)+V^^(l|0)(l-F(O|$)) + ( ,5)mX ^ r ^ ( 1 1 5) 2 9 b (l-$) ra ‘ b . (21) 



5=1 



Finally, using Little (1961) 's theorem the mean packet delay in steady 
state, D, is given by 



D- QymX = { 1-F(0| Q) )S/mXT . 



( 22 ) 



The mean packet delay, D, includes (i) the elapsed time from the arrival of a 
packet until the packet is promoted to the top of the queue, (ii) a number of 
sub (busy)per iods during each of which the packet is not transmitted or is 
involved in a collision, and (iii) the subperiod in which the packet is 
successfully transmitted. 



5. Numeric and Simulation Experiments 
The algorithm to compute D given parameters ( p , m , X , i , 2 ' > was encoded in 
Fortran 77; it first finds 0 and F(0|9) satisfying the condition (14) with 
*®1.E-10, and then compute D using Equations (15) through (22). 
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Computational results are shown in Tables 1 and 2, respectively, for two 
different cases. The two cases have the same values for m (-50) and i' (-3), 
but different values for X and t. The latter parameters were set to 1.6E-4 and 
75, respectively, in Case 1 to represent a high utilization of the channel 
capacity; and to 1.2E-4 and 25, respectively, in Case 2 to represent a low 
utilization. Note that by Lemma 3 an upper bound on the channel capacity 
(i.e., the maximum achievable throughput rate in equilibrium) can be given by 
1/i. The ratio of the total input rate, mX , to 1/i is 62% in Case 1 and 15% in 
Case 2. 



Tables 1 and 2 about here 



In view of the conjecture in Section 4, the range [p^( 1) ,p^ (m) ] of p that 
guarantees stability is provided in each table. Among the steady-state 
measures reported in the tables, the new notation, B + , signifies the set [ b : 

(l\b)> 0). The p-values, at which (l\d)< 0 and B + is nonempty, yield a 

finite mean packet delay, though not satisfying the sufficient condition for 
stability given in Lemma 3. It should be noted, however, that the delay at 
p=*Py(m) is close to the minimum delay: in Case 1 the minimum delay was 128.6 
at the optimal p-value .16, while the delay at p=p^ ( 50 ) - . 08 was 140.4, 
exhibiting a 9% gap from the minimum; in Case 2 the gap was only 2%. It is 
also notable that the expected number of busy users in steady state, or E(|£|) 
(= 7720), was 1.5 and 1.1, respectively, in Cases 1 and 2 when p was set at the 
optimum. This indicates that there is little interference among queues in 
steady state if p is optimized with respect to the mean packet delay, which 
justifies the M/G/l approximation. Indeed, the ratio of collisions to all 
packet transmissions observed in simulation runs was about 5% and 1% in Cases 
1 and 2, respectively, with the optimal p-values. 

The simulation model was written in SLAM II (Pritsker 1984) and run on an 
IBM 3033AP mainframe. In the simulation model, we generated packet arrivals 
for each user so that the interarrival times follow an exponential 
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distribution with intensity A. It is well known that the geometric distribu- 
tion of interarrival times can be closely approximated by the exponential 
distribution when A is small (Feller 1968). This exponential approximation 
tremendously reduced the computing time. Approximately 5,000 packet transmis- 
sions were collected to estimate the mean packet delay for each p-value. The 
first 500 observations were discarded to eliminate the initial transient 
period (see Wilson and Pritsker 1978, Schruben and Goldsman 1985). 

The simulation results are shown in Figure 2 along with the numerical 
approximation results. The solid line and dotted line show approximation 
results for Cases 1 and 2, respectively. Simulation results for Cases 1 and 2 
are plotted using symbols '+' and respectively, and show good agreement 

with approximation results. The delay (versus p) curve is U-shaped and the 
minimum delay is attained in the neighborhood of the maximum p-value that 
leads to equilibrium (viz., near the bottom right corner of the U curve). A 
similar, yet different, pattern of the delay curve was found for slotted ALOHA 
systems by Sidi and Segall (1983). The CSMA/CD achieves close -minimal mean 
packet delays over a much wider range of the p-value than the slotted ALOHA, 
presumably due to the additional medium access control mechanisms — carrier 
sensing and collision detection. 



Figure 2 about here 



6 . Cone lus ion 

A multidimensional queueing process with interactions among individual 
queues arises in many computer and communication systems such as coupled 
processors (Szpankowski 1988), ALOHA satellite communication (Sidi and Segall 
1983) and CSMA/CD local area networks. However, there are no analytic results 
for computing the mean queueing delay in such a system with more than two 
queues . 
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In this paper we have shown that the mean packet delay in CSMA/CD 
networks can be obtained quite accurately using an M/G/l approximation. The 
success of this approximation method is ascribed to the following: (i) The 
message queueing process has the bounded homogeneity property (Lemma 2). With 
this property, after introducing the approximate distribution of the number of 
busy users (Equation (8)), we could simplify the stationary PGF of the packet 
backlog (Equation (7)) into a product form (Equation (9)). (Note that in 
Equation (9), E^_^P(2-q) in Equation (7) is factored out as {F(s | 0 ) -F(0 | 0 ) } . ) 
(ii) The delay analysis incorporates an exact PGF for the random drift of the 
packet backlog over a period between two consecutive embedded Markov epochs 
(Equations (2), (6) and Lemma 1). (iii) Finally, the interference among queues 
is actually minimal in steady state if the control parameter p is optimized. 

A practical and easy way to determine the p-value would be to set it 
equal to p^(m) ; this guarantees stability (by Lemma 3 and Conjecture) and at 
the same time yields a close -minimal mean packet delay in equilibrium (as 
numerically demonstrated in Section 5). According to Conjecture, for given 
system parameters (m, p^(/n) is obtained by solving G^\l\m)= 0 for p. 
Substituting m for b in Equation (A5) in Appendix II, we obtain 

G (1) ( l|ra) ["a ( 1 - p) m +(niAi - l)p(l-p) m *^+Ai' ( 1 - ( 1 - p) m -mp(l- p ) m ~ 1 ) j . (23) 

1-(1 -P) m L J 

Rewriting (23) we see that G^^(l|n 2 )- 0 is equivalent to 

/(<?):- -c 1 q m +c 2 q m ' l +c 3 = 0, (24) 

where q-l-p, c x - c 2 +A(2'-l), c 2 « mA { 2 - 1! ' } - 1 and c 3 = A2'. The equation /(q) has 

the global minimum at q =(m-l) c 2 /mc l . It has exactly two solutions (q x ,q 2 ) 

such that 0< q x < q < q 2 < 1 provided that m\£< 1 (which is a necessary 

condition for stability by Lemma 3) and f(q )< 0. A bisecting search over the 
* 

range (0,q ) of q can quickly find q x , and then we have p^.(ra)= 1-q,. 
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Our approach -- viz., incorporating an exact transient analysis of 
multiple interacting queues into an M/G/l approximation for evaluating average 
delay in steady state -- can be employed to analyze other CSKA/CD protocols. 

It would be interesting to attempt applying this approach to the analysis of 
new variants of CSMA/CD that are designed for high-speed broadcast bus 
networks (e.g., p^-persistent protocol by Mukherjee and Meditch (1988), LCSMA - 
CD by Maxemchuk (1988), and a modified CSMA/CD by Lin and Sousa (1990)). 
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Appendix I. The Conditional Joint PGF of Given 0^*0 

The joint probability of and X(r)« [x^ , iel : 0<x^<r] conditional on Q k 

is given by 



?(R k ~r,X(r)-[x i ,ieI] |Q k )- ^l-(l-p) 



\B k \ + \B' k (r)\ 



HjpM < 



{ || (1-A) r }{|| ) (1-A) W * 1 A(1 -p) r ' W 

^iGl -B k ~B k (r) J ^ie3^( r) w £ x 



r-w 

x r l 



x.-l r-v-x.+l 

A (1-A) 



}• 



(Al) 



In this equation, w represents the slot in which the first packet since the 
beginning of the kth subperiod arrives at user i (who had no waiting packets 
in the beginning of the subperiod). It should be noted that B^(r) is uniquely 
determined given B ^ and X(r) , and B ^ is uniquely determined given Q^. Equation 
(Al) thus shows that the joint distribution of R ^ and X(P^.) depends only on 
B apart from the system parameters p, A, and I. That is, the joint 
distribution is identical for any k and k' with even if , . 

The probability of having a collision in the subperiod beginning at the 
kth epoch, given is 

f |B,|+|B'(r)| 

P(Z k =0|B k ,B£(r))- |l-(l-p) * k - ( | B k | + 1 3 k (r) | ) 

f \B k \ + \B' k (r)U 

|(l-(l-p) K * 

The probability of having a successful transmission in 
at the kth epoch, conditional on k^=r, is simply given 
P(Z^=0 |B^,B^(r) ) . As such, given ^=r, Z ^ depends only 
that B^(R^) is uniquely determined by B , and X(P^) and 
jointly dependent on B^. Therefore, Z^ is transitive!]/ 



P( 



\ D k\ \ D k ^ r ' ‘ “ M / 
1-P) ^ k )/ 



(A2) 



the subperiod beginning 
by P(Z k ^0|B k ,B^(r))- 1- 
on B k and Bj^(r) . Recall 
that R ^ and X(R^) are 
dependent only on B 
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Since X(L^) depends only on Z^, it follows that the random drift depends on 

only through 5^. We thus have proved both Lemmata 1 and 2. 

Let F^+-]_(s iQ/c*®) » w ^ e re s-[s^,iGl], be the conditional PGF of given 

Q^O, and G(s|B):- G k (s|B^-B) be the conditional PGF of given B^-B which is 
independent of k and Q^. (by Lemma 2). Then vre have 



W*IV°>- E 
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Y ) G(s|B) (f^sis^O Vi$B)- \ F ;c (s:s i ,-0,s i -0 Vi$B)l . (A3) 

P( Qjc* 0 ) b-l VBCI:|B|-b 
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Using (Al) and (A2) , function G(s|s) is evaluated as; 
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(A4) 



22 



the 



where II(s)- I | (As ; + l-Ap and $(s)- || (As. + l-Ap are the PGF's of 

iel 1 iel 1 

duration of a successful and an unsuccessful transmission period, 
respectively. A detailed derivation of Equation (A4) can be found in Park and 
Bartoszynski (1990b) . 



Appendix II. Derivatives of the PGF G(s\b) Evaluated at s-l 
Taking the first and second derivatives of G(s|b), which is given in 
Equation (6), with respect to s and evaluating the derivatives at s-1 we 
obtain 



c^dl 



m-b /• ^\ f 

|b)- S tf[{l-(l-p) b+J *}jp(l-A)/(p-A) + (Affli-l)a + 
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(Ami-l)a + Ami' ( 1- (l-p) b+J -a 



A <f> + jp(l-A)V>/(p- A) 
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+ a-a-p) b+j ) 



A (2 n+<y) + 2JpA(l-A)^/(p-A) 



+ j{l-(l-p) fa+j )^j-l + 2jA(l-p)/(p-A) + ( j+1 ) A 2 ( 1 -p) 2 / (p- A ) 2 j | , (A6) 



where a:- (b+J)p(l-p) b+ ^‘ 1 ; 0:- 1- (l-p) b+n (l-A) m * b ' n ; V>:- ) p | ( - 1 )- J ' n /0; 

n= 0 
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n=0 



n=0 
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and 



1k> 

n-0 



(-l)J‘ n (b+n) 2 (l-£) 2 //3 3 . 



Appendix III. Derivation of the PGF r(s|h) and Its First Derivative 
In this Appendix, a, and rfi defined in Appendix II are used again. Based 
on Equation (17) we can write 





s r [i-a-p) b+J - 



a } 



( l_ A ) r ( ra -b-j)(i_ p) rb x 



X 



^ (1-A) 
-w= 1 






(A7) 



Evaluating the summation over w and then the infinite series over r, we obtain 
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Table 1. Steady State Measures for Case 1 where m-50 , i-75 , St ' -3 , A-l . 6E-4 

(p L (l)-l . 6344E-2 , p u (50)-8.4378E-2) 



p 




B + 


E(|5|) 


v (1) a\t) 


F (1) (1 \d) 


r 


s 


D 


1.144E-2 




1 












* 

00 


1.308E-2 




1 


2.2886 


- .1067 


3.9976 


98.32 


427.39 


491.15 


1.634E-2 




- 


2.0754 


- .1370 


2.9448 


91.44 


296.99 


357.26 


1.961E-2 




- 


1.9466 


- . 1598 


2.4314 


86.31 


233.08 


291.19 


2.288E-2 




- 


1.8611 


-.1773 


2.1347 


82.41 


195.91 


252.57 


2.768E-2 




- 


1.7778 


-.1964 


1.8762 


78.13 


163.28 


218.49 


3.902E-2 




- 


1.6721 


- .2242 


1.5871 


71.82 


126.27 


179.65 


5.036E-2 




- 


1.6194 


-.2398 


1.4565 


68.11 


109.17 


161.73 


6.170E-2 




- 


1.5887 


- .2494 


1.3832 


65.63 


99.31 


151.51 


7.594E-2 




- 


1.5655 


- .2570 


1.3273 


63.43 


91.55 


143.64 


8.438E-2 




- 


1.5563 


-.2601 


1.3047 


62.42 


88.29 


140.43 


1.097E-1 


39 


- 50 


1.5405 


- .2655 


1.2611 


60.14 


81.63 


134.19 


1.434E-1 


29 


- 50 


1.5337 


-.2679 


1.2312 


58.00 


76.36 


129.89 


1.603E-1 


26 


- 50 


1.5336 


- .2679 


1.2223 


57.14 


74.49 


128.62 


1.727E-1 


25 


- 50 












* 

oo 



* The system is unstable with V^^(l|5)> 0. 
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Table 2. Steady State Measures for Case 2 where n^50 , i-25 , £ ' -3 , A-l . 2E-4 

(p L (l)-4 . 8303E-2 , Pu (50)-1 . 0722E- 1) 



p 




B + 


E(|B|) V (1) (l\6 ) 


F (1) (l |$) 
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2.415E-3 
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2.4953 
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6 . 5053 
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1065.35 


3.381E-3 
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2.1872 
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Legend: A generation of a new packet 

x transmission attempt sampling a random number greater than or equal to p 
p transmission attempt sampling a random number smaller than p 
Q .total packet backlog in the system at the kth embedded Makov epoch 
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Figure 1 . Channel States and Queueing Process in a Slotted p - Persistent CSMA/CD System 



Figure 2. Comparison between Approximation and Simulation: D versus p 
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